Evolution of twist-shear and dip-shear during X-class flare of 13 
December 2006: Hinode observations 

Sanjay Gosain and P. Venkatakrishnan 
Udaipur Solar Observatory, Physical Research Laboratory, Udaipur, Rajasthan- 313 001, 

o ■ 

(N . INDIA 






6 



> 

(N 

O 

(N 

o 
o 



% 



Received ; accepted 



-2- 



ABSTRACT 



The non-potentiality (NP) of the solar magnetic fields is measured tradi- 
tionally in terms of magnetic shear angle i.e., the angle between observed and 
potential field azimuth. Here, we introduce another measure of shear that has 
not been studied earlier in solar active regions, i.e. the one that is associated 
with the inclination angle of the magnetic field. This form of shear, which we 
call as the "dip-shear", can be calculated by taking the difference between the 
observed and potential field inclination. In this Letter, we study the evolution 
of dip-shear as well as the conventional twist-shear in a (5-sunspot using high- 
resolution vector magnetograms from Hinode space mission. We monitor these 
shears in a penumbral region located close to fiare site during 12 and 13 December 
2006. It is found that: (i) the penumbral area close to the fiaring site shows high 
value of twist-shear and dip-shear as compared to other parts of penumbra, (ii) 
after the fiare the value of dip-shear drops in this region while the twist-shear in 
this region tends to increase after the fiare, (iii) the dip-shear and twist-shear are 
correlated such that pixels with large twist-shear also tend to exhibit large dip- 
shear, and (iv) the correlation between the twist-shear and dip-shear is tighter 
after the fiare. The present study suggests that monitoring twist-shear during 
the fiare alone is not sufficient but we need to monitor it together with dip-shear. 

Subject headings: Sun:sunspots-Sun:flares-Sun:magnetic topology 
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Introduction 



The deviation of an active region (AR) magnetic field from its potential configuration 
can arise due to photospheric footpoint motions and/or during flux emergence. However, 
studies with the high-quality data from space and ground based obser vatories ove r 
the 



2000 



2005 



ast decade suggest that the latter is more dominant mechanism 



Demoulin et al. 



Jing et al. 



2006 



2002 



Falconer et al. 



Schrijver 



2002 



Leka fc Barnes 



2003 



(IWheatland 



Schrijver et al. 



20071 ). The deviation of magnetic field from potential 



configuration is called non-potentiality (NP) of the field. The magnetic energy of an a ctive 



regio n (AR) in excess of the potential magnetic energy is called free magnetic energy ( iLow 



19821 ). The free magnetic energy of an AR, or a portion of it, is released when flare and/or 
C oronal Mass Election s (CMEs) is triggered due to the instability or non-equilibrium 



( jPriest fc Forbes 



2OO2I ). The energy released in the flare is then limited by the amount of 



free energy or magnetic NP of the AR. Therefore, it is important to characterize the NP of 
the AR magnetic field in order to predict the intensity of the flare and/or CME. 

Traditionally, the so-called magnetic shear i.e., the angle between the observed [i/Jo) 
and the potential {ipp) field azimuths, is used to characterize the NP of the active region 
magnetic field and is measured as 



A-i/^ = acos 



B?-BP 



|B?||BP| 

where B° and Bj' are the observed and potential transverse field vectors. The term 
"magnetic shear", used here and in similar studies on solar-magnetism, really refers to 
"observational shear", which is not to be confused with the actual "shear" -a microscopic 
property in fluid dynamics. In what follows, we will use the term magnetic shear for 
traditional reasons and by this term we would always mean the obseryation al shear. Various 



forms of this param eter are used namely, mean shear 



shear (jWang 



19921 ). spatially averaged signed shear (JTiwari et al. 



(IHagvard et al 



1984J ). weighted mean 



20091 ). most probable 
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shear (JMoon et al.l l2003l ) etc. It is well known that the polarity inversion line (PIL) in 
active regions bearirig high ly sheared magnetic fields are the potential sites for flares 



( JHagyard et al. 



1984 . 



19901). These PILs are generally c haracterized by dark fllaments in 



the images taken in chromospheric hydrogen alpha line (jZirin fc Tanaka 



19731). 



It may be noted that all of the shear parameters mentioned above measure the 
twist-shear component of the shear which is measured in the horizontal plane. However, 
one can also measure the shear of a magnetic fleld in the vertical plane. This type of shear 
can be called as the "dip-shear" (we choose the term "dip" because it is synonymous to 
the dip angle measured in geomagnetism), and can be measured by taking the difference 
between the observed (70) and the potential field inclination angle (7^) i.e., A7 = (70 — 7^). 

Physically, the dip-shear can be understood in terms of azimuthal currents, in the 
same way as the twist-shear is understood in terms of axial currents. In our knowledge, 
the dip-shear has not been studied earlier in solar active regions. Henceforth, we will call 
the parameter Aip and A7 and as twist-shear and dip-shear, respectively. The larger the 
value of these angles the larger will be the NP of the AR. It may be noticed that, unlike 
twist-shear, the dip-shear is not affected by the 180 degree azimuth ambiguity, provided 
that the active region is observed close to the disk center. This is because the dip-shear 
depends up on the inclination a ngle of the magnetic field which can be measured without 



ambiguity ( JMetcalf et al. 



20061 ) 



In this letter, we study the evolution of twist-shear and dip-shear in a penumbral region 
located close to the fiaring site in AR NOAA 10930. We use sequence of high-quality vector 
magnetograms observed by the Hinode space mission. This active region was in a (J-sunspot 
configuration which led to a X3.4 fiare and a large CME during 02:20 UT on 13 December 



2006. The fiare was quite powerful and the white light fi are ribbons a. 



lateral motion of the penumbral filaments were observed (JGosain et al 



ong w ith impulsive 



20091) 
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Evolution of the twist-shear and dip-shear together shows interesting patterns which 
can be distinguished in the pre-flare and post-flare stages. In general we find that (a) the 
regions with high twist-shear also exhibit high dip-shear, (b) the penumbral region close 
to the flare site shows high twist-shear and dip-shear, and (c) twist-shear and dip-shear 
studied together can be used to study flare related changes in the active regions. 

The paper is organized as follows. The observational data and the methods of analysis 
are described in section 2. The results are presented in section 3 and the discussions and 
conclusions are made in section 4. 

2. Observational Data and Methods of Analysis 
2.1. Hinode observations 



A sunspot with 6 configuration was observed in AR NOAA 10930 during 12 - 



13 December 2006 



Ichimoto et al. 



3y the Spectro-Polarimeter (SP) instr ument ( 



2008) with So: 



ar Op tical Telescope (SOT) ( iTsuneta et al 



Lites et al. 



2007 



20081 ) onboard 



Hinode satellite (iKosugi et al.ll2007l ). The SP obtains the Stokes profiles, simultaneously 
in Fe 1 6301.5 A and 6302.5 A line pair. The spectro-polarimetric maps of the active 
region are made by scanning the slit across the field-of-view. This takes about an hour 
to complete one scan. We choose a sequence of six SP scans from 12 December 2006 
03:50 UT to 13 December 2006 16:21 UT when the sunspot was located close to the disk 
center with heliocentric distance (//) of 0.99 and 0.97, respectively. The scans were taken 
in "Fast Map" observing mode with following characteristics: (i) Field-of-view (FOV) 295 
X 162 arc-sec, (ii) integration time of 1.8 seconds and (iii) pixel-width across and along 
the slit of 0.32 and 0.29 arc-sec, respectively. Th e Stokes profiles were then fit ted to an 



analytic solution of Unno-Rachkovsky equations (jUnno 



1956 



Rachkovsky 



19731) under the 
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assumptions of Local Thermo dynamic Equili 



jrium (LTE) and Milne-Eddington model 



atmosphere (ILandqlfi fc Landi 



code called HeliX (JLagg et al 



Deg^ 



'Innocent! 



19821 ) with a non-linear least square fitting 



2004J ). The physical parameters of the model atmosphere 
retrieved after inversion are the magnetic field strength, its inclination and azimuth, the 
line-of-sight velocity, the Doppler width, the damping constant, the ratio of the line center 
to the continuum opacity, the slope of the source function and the source function at r 
= 0. We fit a single component model atmosphere along with a stra y light component. 



The inversion code HeliX is based upon a reliable genetic algorithm ( ICharbonneau 



19951). 



This algorithm, although slower, is more robust than the classical Levenberg-Marquardt 
algorithm i n the sense that the global minimum of the merit function is reached with higher 



reliability (JLagg et al. 



2004h . 



The vector magnetograms obtained after inversi on, were first solved for 180 degree 



azimuth ambiguity by using the acute angle method (JHarvey 



19691 ) and then transformed 



from the observed fra me (image plane) to loca l solar frame (heliographic plane) using the 



procedure described in 



Venkatakrishnan et al. 



( 1l988l ). The potential field wa s computed 



from the line-of-s ight field component by using the Fourier transform method (jAlissandrakis 



1981 



Gary 



19891 ). The IDL routine used for potential field computation is fff .pro which 



is available in the NLFFF (Non-linear force free field) package of the SolarSoft library. The 
continuum intensity images of the sunspot, corresponding to the sequence of scans used, are 
shown in figure 1, with the transverse field vectors overlaid on it. The two magnetograms 
were aligned using the cross-correlation technique applied to the continuum image of the 
sunspot. A black rectangle is overlaid on these images to show the location of the region 
where we study the evolution of twist-shear and dip-shear. The location of this black 
rectangle is chosen with the help of a co-aligned G-band filtergram observed from Hinode 
Filtergraph (FG) instrument during fiare. This G-band image is shown in the left panel 
of figure 2. The fiare ribbon is marked by '+' symbols and the black rectangle of figure 
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1 is shown here also. The flare ribbons sweep across the rectangular box during 02:20 to 
02:26 UT. This indicates that the rectangle is chosen such that it samples the penumbra 
which is very close to the flaring site. The right panel of flgure 2 shows the longitudinal 
magnetogram in order to indicate the location of rectangle (shown here with white color) 
with respect to the PIL. The maps of dip-shear A7 and twist-shear Aip for the sequence 
of vector magnetograms are shown in flgure 3 and 4 respectively. In flgure 5 we show the 
distribution of dip-shear A7 and twist-shear Aifj inside the black rectangle and its evolution 
with time. 

3. Results 
3.1. Distribution and Evolution of Non-Potentiality in NOAA 10930 

3.1.1. Dip- Shear 

The maps of dip-shear A7 for the entire sequence of vector magnetograms covering 
the pre-flare (panels (a)-(c)) and post-flare (panels (d)-(f)) phases are shown in flgure 3. 
The value of fleld inclination 7 is measured with respect to local solar vertical direction and 
ranges from to 180 degrees. For purely vertical positive (negative) polarity fleld the value 
of 7 corresponds to degrees (180 degrees). The black rectangle in flgure 3 corresponds to 
negative polarity fleld. Therefore, the positive value of dip-shear A7 inside this rectangle 
means that the observed fleld is more vertical than potential fleld. The magnitude of 
dip-shear A7 can be judged with the aid of colorbar at the bottom of flgure 3. 

It may be noticed that: 
(i) the value of dip-shear A7 is large inside the rectangle as compared to other penumbral 
locations, 
(ii) in the pre-flare (panels (a)-(c)) phase the dip-shear A7 consistently has a large 



magnitude which decreases in the post-flare phase (panels (d)-(f)). 



3.1.2. Twist- Shear 



The field azimut h has been so lved for 180 degree azimuth ambiguity by using the 



acute angle method (jHarvey 



19691 ) a nd the projection effects hav e been removed by the 



application of vector transformation (jVenkatakrishnan et al. 



19881 ). This azimuth ambiguity 



resolution method works well in the regions where the angle Aifj is less or greater than 90°. 
For regions where Aifj reaches value close to 90° such as along parts of polarity inversion 
line (PIL) in fiaring active regions, the accuracy of the method is poor. This is why we 
choose a rectangular box for studying the evolution of twist-shear and dip-shear sample 
the penumbra close to fiaring site and at the same time stay away from the PIL where the 
acute angle method may have problems in resolving the azimuth ambiguity. 

The maps of twist-shear Aip for the entire sequence of vector magnetograms covering 
the pre-fiare (panels (a)-(c)) and post-fiare (panels (d)-(f)) phases are shown in figure 4. 
The value of field azimuth ■0 is measured with respect to the positive X-axis and is positive 
in the anti-clockwise direction. The magnitude of twist-shear A-^ can be judged with the 
aid of colorbar at the bottom of figure 4. However, for the present study the sign of shear 
angle is not important, so we shall focus on its magnitude. It may be noticed that the 
value of twist-shear Aifj is large inside the rectangle and adjacent PIL as compared to other 
penumbral locations. However, the fiare related changes are not so discernable to eye as 
compared to dip-shear in figure 3. 
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3.1.3. Evolution of Dip- Shear and Twist-Shear 

In figure 5 we show the scatter between the dip-shear A7 and twist-shear A^ for the 
pixels within the black rectangle shown in previous figures. The panels (a) to (c) correspond 
to pre-fiare while the panels i.e., (d) to (f) correspond to post-fiare phase. 

It may be noticed that : 
(i) The distribution of dip-shear A7 and twist-shear Aip in panels (a)-(c) is different from 
the distribution in panels (d)-(f). 

(ii) The dip-shear increases before the flare (panels (a)-(c)) but twist-shear tends to decrease 
at the same time. 

(iii) The dip-shear and twist-shear are in general correlated, i.e. the pixels with large 
dip-shear also have large twist-shear. 

(iv) The most important change can be noticed after the flare, i.e., between panels (c) 
and (d). After the flare, (panel (d)) the dip-shear decreases significantly while twist-shear 
increases. However, now both shear components show less dispersion i.e., follow a tight 
correlation. 

(v) In panels (d)-(f) the two shears maintain smaller dispersion but dip-shear starts to 
increase once again. This increase suggests that the non-potentiality was building up again 
in the active region. It may be noted that the flaring activity continued in this region on 
the next day i.e., 14 December 2006 also, with another X-class flare occurring at about 
22:00 UT. 

4. Discussion and Conclusion 



HudsonI (I2OOOI ) conjectured that the free-energy is stored in non-potential magnetic 



loops that are stretched upwards and the free-energy release during the flare must be 
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accompanied by sudden shrinkage or implosion in t 



le field. Also, it is p redicted that after 



the flare the field should become more horizontal (JHudson et al. 



20081 ). Using coronal 



images durin g flares, there ar e observational reports about detection of the loop contraction 



during flares (JLiu et al. 



20091 ). 



Further, it was shown by i Venkat akrishnaru (jl990l ) that in force-free fields a high 



non-potentiality implies weaker magnetic tension, which in turn implies a larger vertical 
extension of the field due to lower magnetic pressure gradient. Conversely, the release of 
free magnetic energy during flare implies a loss of magnet ic non-potentiality lea ding to a 



decrease in the vertical extension of the field or shrinkage (JForbes fc Acton 



19961). 



The non-linear force 
region by ISchrijver et al. 



ree f ield (NLFFF) extrapolations of the NOAA 10930 active 



( 120081 ) show the non- potentiality of this activ e reg ion in the forrn 



Venkat akrishnani (jl990l ) and 



Hudsonl teoOOh 



of a twist flux rope structure. As suggested by 
such a structure will have larger vertical extension in pre-flare as compared to the post-flare 
configuration. The closer the post-flare field approaches to the potential field configuration 
the smaller is the value of inclination difference A7 expected. This may give an explanation 
for the decrease of dip-shear A7 after the flare. 

However, in contrast, the increase in the twist-shear ^ip after the flare also needs an 
explanation. The opposite behaviour of twist-shear and dip-shear in relation to the flare 
can be understood in the following way. The twist-shear is dependent on sub-photospheric 
/ photospheric forces, so the twist-shear will continue to increase independent of coronal 
processes like flare. However, the plasma j3 decreases rapidly above the photosphere and 
thus there is no non-magnetic force or shear that is strong enough to change the inclination 
of the field lines. Hence inclination will be more responsive to coronal processes. This may 
explain why inclination became more potential after the flare. Hence, dip-shear could be a 
better diagnostic of NP above the photosphere. 
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In summary, we studied the evolution of twist-shear and dip-shear in a flaring 
(J-sunspot using a sequence of high-quality vector magnetograms spanning the pre-flare 
and post-flare phases and found that : (i) the penumbra located close to the flaring site 
has high twist-shear and dip-shear as compared to other parts of the penumbra, (ii) the 



Jing et al. 



(120081) also, (iii) 



twist-shear increases after the flare which was earlier reported by 
the dip-shear however shows a decrease after the flare, (iv) the twist-shear and dip-shear 
are correlated, i.e. pixels with high twist-shear exhibit high dip-shear, and this correlation 
is much tighter after the flare, and (v) distribution of twist-shear and dip-shear and its 
evolution (in figure 5) clearly shows different patterns before and after the flare. 

This type of behaviour in the twist-shear and dip-shear parameters will need to be 
evaluated further in more flares before it can be understood physically. We plan to carry 
out more extensive study of the dip-shear and twist-shear in existing Hinode datasets. 
However, a high-cadence study of these shear parameters would be possible only with 



the upcoming observations fror n Helioseismic and Magnetic 



Dynamics Observatory (SDO) (iScherrer fc SDO/HMI Team 



mage r (HMI) onboard Solar 



2OO2I ). The present study is 



important in the sense that it points the way to a vector-field follow-up to the results of 



Sudol fc Harveyl (120051 ) . which established the line-of-sight field changes during powerful 



flares. 

In the context of present study, one should bear in mind that the vector magnetograms 
derived from the Hinode SOT/SP scans, although polarimetrically very precise are very 
noisy geometrically. An unwanted consequence of the geometric noise could be that the 
flows, specially on long time scales, would tend to create an appearance of non-potentiality, 
even if there was none. This is an important issue which needs to be addressed sooner 
than later, considering the widespread use of SOT/SP magnetic maps as the "vector 
magnetograms" . We plan to carry out a detailed study of this effect using simultaneously 
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observed Spectro-Polarimeter (SP) scan from Hinode SOT and vector-magnetograms from 
HMI onboard SDO. 

We thank the anonymous referee for his/her valuable comments and suggestions, 
specially for pointing out the geometric noise in SOT/SP scans and its side effects. We 
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the post-launch operation is provided by JAXA and NAOJ (Japan), STFC (U.K.), NASA, 
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Fig. 1. — The panels (from top to bottom) show continuum intensity map of the 5-sunspot 
in NOAA 10930 during the times mentioned at the top. The transverse magnetic field vectors 
are shown by arrows overlaid upon these maps. The black rectangle, shown in all panels, is 
the region where we monitor the evolution of twist-shear and dip-shear. 
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Fig. 2. — The left panel shows the G-band filtergram of the 5-sunspot in NOAA AR 10930 
during 12 December 2006 02:30 UT, with the location of the flare ribbon marked by '+' 
symbols. The flare ribbons sweep across the rectangular box during 02:20 to 02:26 UT. The 
right panel shows the map of the longitudinal field component for this sunspot. The black 
(white) rectangle in left (right) panel marks the region where we monitor the evolution of 
twist-shear and dip-shear. 
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Fig. 3. — The panels (from top to bottom) show the maps of dip-shear, A7 = Jo ~ Ip for 
NOAA AR 10930 at different times. The black rectangle marks the region where we monitor 
the evolution of twist-shear and dip-shear. 
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Fig. 4. — The panels (from top to bottom) show the maps of twist-shear, Aip = ipo ~ i^p for 
NOAA AR 10930 at different times. The black rectangle marks the region where we monitor 
the evolution of twist-shear and dip-shear. 
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Fig. 5. — The panels (from top to bottom) show the evolution of twist-shear A-^ and 
dip-shear A7 inside the region marked by rectangle in previous figures. 



